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1. Introduction 

In system with strong electron-phonon (e-pli) interaction, the carriers lose mobil- 
ity, ultimately acquiring polaronic character. A polaron is a state in which the phonon 
and electron degrees of freedom are strongly entangled, and the presence of an electron 
is associated to a finite lattice distortion, which in turn binds the electron leading to 
the so-called self-trapping effect. Polarons also tend to create bound pairs, called bipo- 
larons of course the presence of Coulomb repulsion destabilize bi-polarons in favor of a 
pure polaronic state at finite densities ^E]- Typical signatures of polarons are seen in 
multi-peaked photoemission spectra[3] and transport measurements, where an activated 
behavior with a characteristic energy given by the polaronic binding energy is observed. 
The polaronic peak found in the mid infrared measurements of optical conductivity 0] 
may also not only detect the polaronic binding energy |^ but also other subtle polaronic 
transitions at very low energy p]. Another less classical indication of polaronic forma- 
tion comes from the analysis of lattice displacements associated to the excess charge as 
obtained by the distribution of distances between atoms|3j. 

A joint analysis of both spectral [E] and local lattice distortions can disentangle various 
kind of behavior in polaronic systems. In fact we notice that neglecting the repulsive 
interaction a gapped pair state can be formed even without a significant associated 
polarization. In this case bipolarons are expected to be a relatively mobile particle. 

The aim of this work is to provide a thorough analysis of polaronic spectral properties 
in both single and multi-polaron cases. The backbone of our presentation is the DMFT 
which is introduced and briefly discussed in general in section [2l 

In the single polaron case we review the exact solution of the Holstein model [3 |H] 
(section ^ and we present some new results for the Holstein t — J model comparing 
our results with those of Mishchenko "H* . At large polaronic densities we use both Exact 
Diagonalization (ED) and Quantum Monte Carlo (QMC) techniques to solve the DMFT 
equations for the Holstein model respectively at the T = (sec. |5J| and T > (section 
l6|l . In this case we compare spinless and spinful fermions cases and we discuss in detail 
the role of the adiabatic ratio. In this way the properties of a pure polaronic state can 
be disentangled from those of bipolaronic state. We compare also numerical solutions 
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with analytic approximate schemes based on Born Oppenheimer (|5'l|l and Lang-Firsov 
canonical transformation l|5'2|l respectively in adiabatic and antiadiabatic regime. 

2. The DMFT method 

2'1. Introduction. ~ Dynamical Mean Field Theory is a non-perturbative technique 
originally developed as the exact solution of a interacting electron problem on an infinite 
dimensional lattice ^3 m . A comprehensive review can be found in ref . ^| now let us 
sketch some key point to understand the developments presented in following sections. 
Let a consider a general tight-binding problem on a lattice with coordination number z 



(1) 



where Ci are fermionic annihilation operator acting on site i (spin index omitted for 
simplicity) 1^ is a local (on-site) potential. The scaling of hopping t is such as the limit 
z — > oo give a non trivial result. Mean field theory turns out to be exact in infinite 
dimensions i.e. when the number of nearest neighbor diverges, therefore we can replace 
the effect of hopping on neighboring sites as cavity field rj 



with the sum running on nearest neighbor of site i. In terms of the internal fields r], 
Hamiltonian can be written formally as a sum of single site operators as 



(2) 



^ = - - E ^l^' - E + E ^[^- 4] 



For fermions rj obeys anticommutation relations 



= tS^ 



On a more formal ground it can be demonstrated that the cavity field is a Gaussian 
Grassmann field which is therefore determined solely by its correlation function jl2j. 

2'2. Single impurity action. - The previous arguments can be more formally developed 
using a path- integral formalism. In analogy with classical mean- field T^, the fermions of 
all sites but one (namely 0) are integrated out leading to a single-site partition function 



(3) 



where now ■0 and ■0^ are Grassmann anticommuting eigenvalues of the creation/destruction 
operators of site 0^3- The action S is given by 



(4) S = - dT dT'^^{T)go\T-T')^{T')+ rfTF[0(r)>^(T)] 

Jo Jo Jo 
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where the correlator Qq^{t — t') is 

(5) Q^\r t') = dr5{T t') + {m{r)m{r')) 

which is independent on i due to translation invariance. The action depends para- 
metrically on the environment through the correlator ^o- Q is indeed the action of a 
single impurity embedded in a medium whose properties are related to the original lat- 
tice via the self-consistency equation (0). To be more concrete let us consider an infinite 
coordination Bethe lattice of half-bandwidth D, for which the Green's function of 77 is 
proportional to the local Green's function Gjj- ^21 

j 

then Eq. © reads Gq^t) ^ -dr + {D'^/^)G{t), or in frequency domain 

(6) GQ^ioJn) - ioJn + {D^/^)G{iuJn). 

Eqs © and Q), together with the self-consistency condition form a closed set of 
mean-field equations. 

2'3. Single impurity Hamiltonian. - A Hamiltonian formalism can be also developed 
when suitably defined fermion fields are introduced to get required gaussian cavity field. 
Eq. can be obtained by integrating out auxiliary fermionic fields Ck of an Anderson 
Impurity Hamiltonian (AIM) 

(7) Haim = Mf^ck + cfc/t) + ^ Ekc'^ck + V[f, f^] 

k k 

levels Ek and hybridization constants Vk must be chosen to give the appropriate cavity 
field. In the Bethe lattice case 

where in the l.h.s. we read the local lattice propagator. Eq. (|SJ| becomes the self- 
consistency condition which determines the appropriate AIM parameters Ek and Vk- 

3. Holstein model in infinite dimensions 

The Holstein molecular crystal model is the paradigmatic model for small polarons. 
Its Hamiltonian reads 



H = -tY^ {cl^CJ.a + h.c.)-gY{n,^a-^)ia^ + al) 



(9) + ojo^ ajoi 
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where Ci^a (cj^,) and Oi (a|) are, respectively, destruction (creation) operators for fermions 

and for local vibrations of frequency ujq on site i, rii^cr ~ ^Ci^a the electron density per 
spin, t is the hopping amplitude, g is an electron phonon coupling. In the half-filled case 
we always fix the chemical potential to the particle-hole symmetric value, which fixes the 
density per spin to n = 1/2. In the spinless case there is no sum on a. We choose as 
parameter of the model the c-ph coupling constant A = 2g^ /ujqD where D is the half- 
bandwidth of our infinite-coordination Bethe lattice, and the adiabatic ratio 7 = loq/ D. 
In the single electron case the spin index is unessential, moreover the mean density is 
zero and the el-ph interaction in Eq. (jSJ is replaced by — <?X]i + ^\)- The partition 
function ^ is now defined as 



(10) I Vx{t) I r»V^2'V^(T)exp(-S') 



where using units where the spring constant K = Mujq ~ 1 x{t) = y^ujf)/2{a{T) + a'^{T)). 

The single impurity action S associated to the lattice Hamiltonian in an infinite 
coordination lattice reads this case: 

(11) S = - f dr r dr'J2i:l{r)g^\r-r')Mr') + 

Jo Jo 



2 Jo V ^0 



(12) -I- - / dri —::^+x\t) 



fl3 

(13) -I- VA / dTx{T) (n{T) ~ 1) 



where ^(t) = J2a'^ti'^)'^<^i'^)- 

31. DMFT-QMC method. - An efRcient method to solve Eqs. H11I12I13|I at finite 
temperature is the QMC method (in the Blankenbecler-Scalapino-Sugar approach|14|. 
This method works well for not too low temperatures and naturally yields electronic 
and bosonic correlation functions, as well as the probability distributions associated to 
the phonon fields. The method is not affected by the negative sign-problem, and its 
main limitation comes in the adiabatic regime (7 <C 1) where the phonon becomes heavy 
making more difficult to sample correctly the available phase space. In the ESS scheme 
the fermions are integrated out , and the phonons coordinates x{t) are discretized in 
L imaginary-time slices of width At — f3/L and then sampled by QMC. L has to be 
chosen large enough to reduce as much as possible At, which controls the the Trotter 
discretization error. To keep At less than 1/8 we used 32 slices except for the lowest 
temperature (/3 = 8) for which we have used L = 64, 

3'2. DMFT-ED method. - The Anderson Impurity Model for the Holstein model reads 

HaIM = = - X! ^fc(4,a/'^ + ^■<^) + X! ^kcl,aCk,a " 
k.a k,(T 

(14) - gY^(flf,-^{a + a^)+u,a^a. 
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We solve (I14|l by means of ED by truncating the sums in the first two terms of Eq. 
(|14|) to a small number of terms Nb, so that the Hilbert space is small enough to use, 
e.g., the Lanczos algorithm to compute the T — Green's function. For the case of 
phonon degrees of freedom we consider here, also the infinite phonon space has to be 
truncated allowing for a maximum number of excited phonons Nph- In all the calculations 
presented here the convergence of both truncations have been checked. The value of 
Nph has to be chosen with special care in the adiabatic regime and in strong coupling, 
where phonon excitations are energetically convenient. As far as the discretization of the 
bath is concerned, the convergence of thermodynamic averages and Matsubara frequency 
properties is exponentially fast and iVf, ~ 8 — 9 is enough to obtain converged results. 
The method also offers the advantage of a direct evaluation of real-frequency spectral 
properties such as the electron and phonon spectral functions. The main limitation is 
that these quantities reflect the discrete nature of our system. In practice, the spectra 
are formed by collections of 5- functions, which limits our frequency resolution, and makes 
the method better suited to gain knowledge on the main features of the spectra, rather 
than on the fine details. 

3'3. Quantities of interest. - We mainly characterized electronic and phononic prop- 
erties by considering respectively the electron density of states p{uj) ~ —^G{u!) and 
the phonon probability distribution function (PDF) P(x) = {4>o\x) {x\(j)o) ^ where is 
a phonon coordinate eigenstate, and \(j)o) is the ground state vector. At finite electron 
density a lattice polarization reflects in the presence of two peaks in P{x), corresponding 
to opposite polarization of occupied and unoccupied sites (bimodal behavior) . In the 
single electron case the polarization of a lattice site where one electron sits is a marker 
of a polaronic crossover. Also in this case a bimodal behavior can be observed in the 
adiabatic regime but generally speaking we have a definite polarization when phonon 
fluctuations are less than the average polarization due to the presence of the electron. 
In this way a qualitative difference is identified between the polarized and unpolarized 
regimes, which allows for an unambiguous way to draw a crossover line, as opposed to 
estimates based on smoothly varying functions as average lattice fluctuations or electron 
kinetic energy. 

A Metal to Insulator Transition (MIT) can be probed by the low energy behavior 
of p. This is a consequence of the infinite dimensional limit where the self-energy is 
momentum-independent. Thus the vanishing of the low-energy quasi-particle spectral 
weight coincides with the divergence of the effective mass, which determines the MIT. 



4. Single electron in Holstein and Holstein t-J models 

4T. T ~ continued fraction for a single polaron. - The single electron case in both 
Holstein, t — J and t — J-Holstein case can be solved semi-analytically taking advantage 
of peculiar features of the zero density case. We briefly describe here the formalism at 
T = 0. Generalization to thermalized lattice can be found in We use here the AIM 
formalism of eq. (|14|l . For a single electron the Green's function is purely retarded, then 
the retarded impurity propagator can be defined as 

(15) G(c^) = (01/ „ f^\0) 

u} + Id — H 



which has the correct prescription 6 > for convergence of the time integrals. The 
vacuum energy is defined here to be zero. To proceed further one needs to introduce the 
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generalized matrix elements 



(16) = (0|4=/ , .\ ^ /^^|0) 



UJ + iS - H' 

so that the element Go,q will be the Green function. 

Let us separate the impurity Hamiltonian of eq. (|14|l into Hq and Hj, Hj is the local 
interaction term and Hq the remainder. A useful operator identity for the resolvent is 

1111 
(17) Tr + TtHi- 



z-H z-Ho z-Hq z-H 

The diagonal matrix element of this operator on the impurity zero phonon state /t|0) is 
the Green's function of eq. H15|) . 

In the subspace of zero electron p-phonon states \0,p) = (at)P/-^/pI|0) one can write 

(18) ff/ = ^/t|0,p)(0,p|/(a + at) 



leading to the recursion formula for the Gn,m^s 

(19) Gji^rii — Gon^n^ni S ^ ^ GQn-^n,pGp^rn, 

P 

where Gon = Gq{lo — nuo) is the diagonal element of the free resolvent and Xn^p are the 
phonon displacement matrix elements Xn,p ~ \/p + ISn.p+i + ^/pSn.p^i- 

One immediately recognizes that, due to the particular form of the matrix X, G^^ is 
a tridiagonal matrix, so that the solution of the problem is reduced to the inversion of a 
matrix in arbitrary dimensions. Following the lines given in Ref.^lB) one can express the 
diagonal element of the G matrix in terms of the diagonal and non-diagonal elements of 
G~^. The local propagator (the 0,0 element of G) is obtained in terms of a Continued 
Fraction Expansion (CFE), as a functional of the "bare" propagator Gq: 



(20) G{lu) 



1 



Cq ^(W - L^o) - 



Go\uJ-2u;o) 



33^ 



Go (uj - 3^0 ) 



Due to the impurity analogy, this is also the local propagator of the original lattice 
problem, provided that self-consistency condition (eq. ® for real frequencies) is fulfilled. 
As a consequence the lattice self-energy E is immediately obtained by G = 1/{Gq^ — E). 
An example of the spectral function that can be obtained with this formalism is presented 
also in this volume (see |E| fig- 4). 
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4'2. Holstein t — J model in infinite dimensions.. - To derive the Holstein-t — J 
Hamiltonian for a single hole we sketch the treatment of Ref. ^|. First the t ~ J 
Hamiltonian is transformed by a canonical transformation into a ferromagnetic one. Then 
we introduce fermionic h hole and bosonic h spin-defect on the antiferromagnetic ground 
state operators. As a further step an Holstein-Primakoff transformations is performed 
which introduces spin waves |19| , then we adopt the linear spin wave scheme and we use 
explicitly the infinite coordination limit to get: 

H = ^ {h]hiaj + h.c.^ ~ .^X! ~ i^l^i) (oi + 4) + X! '^1"' 



(21) +i^E[^'^^ + %]+iE'*l^' 



J 



The first term of Eq. (|21() describe the kinetic hopping of one hole on the antiferromag- 
netic background, which is accompanied by the creation (destruction) of a spin defect 
which breaks (restores) 2z magnetic bonds with individual energy J/4z. In addition we 
have the usual local e-ph interaction which couples the hole density to the local phonon. 
The last term in Eq. I|21|l can be absorbed in the definition of the hole chemical potential 
which, for the single hole case here considered, has to be set at the bottom of the hole 
band. The single-hole Green's function is given by the resolvent of Eq. H15() in which the 
spinless fcrmions / are replaced now by the hole operators h. Following the same path 
of the previous calculation and using the results of Ref. [23 for the t~ J model, we can 
write the local hole propagator G{u!) as function of the sum of a hopping and a phonon 
contribution to the self-energy [H^: 

(22) G{uj) = 



— 5]hop(w) — Eel_ph(w) ' 

where the hopping contribution describes the dynamics of the hole through the antifer- 
romagnetic background, 

(23) Ehop(w) = jG(t^- J/2), 

while the e-ph self-energy takes into account the on site multiple scattering with phonons 
and it is formally the same as in Eq. (|2()|l after identifying 

(24) Go-i(a;)=c.-^G(c.- J/2). 

Both Ehop(w) and Eei_ph(w) are expressed as functionals of the total Green's function 
G{u}) leading to a self-consistent interplay between the spin and the e-ph interaction. Eqs. 
(|22|) , (|23|l and H24|) represent a closed self-consistent system which we can be numerically 
solved by iterations to obtain the explicit exact expression of the local Green's function 
G(tf) The formal scheme looks quite similar to the single electron solution of the 

Holstein model However, due to the antiferromagnetic background, the physical 
interpretation is quite different. 

Due to the orthogonality of the initial and final antiferromagnetic background, the 
non-local component of the Green's function in the Holstein-t— J model, even for J = 0, is 
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strictly zero Gij{Lo) — G{Ld)5ij,^^ whereas for the pure Holstein model Gi^j{uj) is finite 
and provides informations about the non-local dynamics: G(k, cj) = 1/1'-^ ~ Ck ~ ^(i^)]- 
In addition, the magnetic ordering has important consequences also on the local 
Green's function Gu^lu). In the pure Holstein model for instance Gii{uj) takes into 
account any generic dynamics which occurs back and forth a given site whereas in the 
Holstein-t-J model the electron must follow a retraceable path in order to restore the 
antiferromagnetic background. [2DI A Bethe-like dynamics is thus enforced by the mag- 
netic ordering regardless the actual real space lattice. The object made up by the hole 
plus the local modification of the spin configuration due to the presence of the hole is 
the "spin polaron" . 

4'3. Results. - The local physical properties of one hole in the infinite dimension 
Holstein-i — J model have been extensively investigated in Ref. [21 ■ In Fig. ^ we report 
the evolution of the local spectral function p[ijj) = — (l/7r)ImG(a;) as function of the 
spin-exchange J /t and of the e-ph interaction A for uja/t = 0.1. In the limits J = and 




Fig. 1. - Hole spectral function in the t — J-Holstein model for 7 = 0.1 as function of J/t and A. 
Left panel: J/t = 0.4 and (from bottom to top) A = 0.0,0.05,0.15,0.25,0.35,0.45,0.55. Right 
panel: A = 0.25 and (from bottom to top) J/t = 0.0, 0.1, 0.5, 1.0, 1.5, 2.0. 

A = we recover respectively the results of Refs. (Sj and |2Qj. They are qualitatively 
different: for J/t = and finite A the spectral density is the same as the pure Holstein 
model (see e.g. fig. 5 ^21 of the present volume), whereas for A = and finite J/t the 
spectrum is described by magnetic peaks spaced as {J/t)'^/'^ (for small J/t). Switching 
on at the same time both the magnetic interaction and the e-ph interaction gives rise 
to the interesting interplay between these degrees of freedom. This can be shown for 
instance in the evolution of the spectral function as function of A: increasing the e-ph 
interaction not only gives rise to additional phonon peaks which superimpose on the 
magnetic one, but it also enhances the nature of the magnetic polaron, from a large one 
(corresponding to ( J/t)^/^ spaced peaks) to a more localized one (corresponding to {J/t) 
spaced peaks). A similar behavior appears as function of J/t: here for large J/t the gross 
structure of the spectral function is described by equally (J) spaced magnetic structures 
with fine features determined by the e-ph interaction. Note that this latter is represented 
by phononic peaks spaced by as in the antiadiabatic limit, although 7 here is just 
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7 = 0.1, since the antiadiabatic regime is intrinsically enforced by the reduction of the 
kinetic energy due to the magnetic polaron trapping. In any case within DMFT we loose 
hole band dispersion. As a consequence we were not able to reproduce even qualitatively 
the interchange of magnetic and polaronic peaks observed by increasing el-ph coupling 
at finite dimensionality However within our localized solution we recover many of 
the qualitative features of both magnetic and lattice polaron crossovers QHI ■ 

Another interesting quantity which points out the interplay between magnetic and 
e-ph interaction is the PDF P{x). Notice due to localization of DMFT solution P{x) 
gives the lattice distortion associated to the localization center. In the left panel of Fig. 
|2we show the P{x) for J/t = 0.4 as function of A. Here the lattice polaron formation 




-2-1 012345 -1 1 2345 

X X 

Fig. 2. - Phonon PDF in the t — J-Holstein model for coo/t = 0.5 Left panel: J/t = 0.4 and 
(from left to right) A = 0.0,0.05,0.25,0.5,0.75,1.0,1.25. Right: A = 0.5 from left to right 
J/t = 2.0, 1.0,0.40.05. 

is characterized by the value Xpoi at which broadening of the PDF is maximum. For 
larger A the PDF recovers a gaussian form due to the lattice fluctuations around the new 
minima at finite distortion. The magnetic polaron formation is also pointed out by the 
analysis of the P{x) (right panel of Fig. [5J: the magnetic trapping favors the lattice one, 
further reducing the anomalous lattice fluctuations towards the gaussian ones. 

5. Half- filled Holstein model: spinless vs spinful cases at T = 

5'1. Adiabatic regime. - We briefly describe, as a starting point of the Born-Oppenheimer 
(BO) procedure, the adiabatic limit in which 7 — > keeping A fixed. This limit has been 
thoroughly studied in Ref. \T5\ which we here briefly resume. When uq and A 
finite the kinetic term l|12|l forces the phonon path x{t) to be r-independent, x{t) = x. 

Phonons becomes classical and the interaction term reads — \/Ax dT{n(T) — 1). The 
gaussian integrals in (|10|l with the action given by (|11I12I13II can be computed analyti- 
cally leading to 



(25) Z= dxexp{-pV{x)) 
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where the adiabatic potential V{x) is |21[I10[|^ 

^''^ ""^'^^r ^^^^^^ j 

where s is spin degeneracy. In formula Eq. (|26|l we have found useful to separate the 
contribution in absence of hybridization ( first line of Eq. p6|) ) from a remainder (last 
line). Through the adiabatic potential we compute the phonon PDF as 

(27) P{x) = ^"P(-f 

Taking advantage of the Gaussian nature of the fermions we get for the local propa- 
gator 

(28) G(w) = / dxP(x 



which defines the self consistency condition through Eq. ©. 

The self-consistency condition © together with Eqs. H28I27I26|I . completely solves 
the problem. Notice also the correspondence between the spinless and and the spinful 
case upon rescaling A to A/2 in the latter case IF. For the Bethe lattice (at zero 
temperature) it can be shown that the potential (|26|l becomes double welled above a 
critical value Apo/ — 37r/(8s). A MIT occurs at a larger couphng \mit — 1.328/s[l5|. 

The BO procedure goes on by quantizing the adiabatic potential after adding the 
phonon kinetic energy contribution. Introducing the scaled variable u = gxj ^ the BO 
phononic Hamiltonian reads 

(29) ij^o^__L|L+,y(u). 

and V(u) is given by Eq. (|26|l . Notice that the spinful BO Hamiltonian maps onto twice 
the spinless one upon rescaling 

A/2 ^ A 

(30) 27 ^ 7. 

While phonon properties are immediately obtained at this stage from the solution of 
the one-dimensional anharmonic system of Hamiltonian Eq. H29f) . electronic properties 
must account non-trivially for the tunneling of phonon coordinates. 

The simplest way to describe electrons coupled to a tunneling system is to map it 
onto a two level system. In our model this can be accomplished by changing the basis 
(operators a) from that of the harmonic oscillator to the that defined by the solution of 
(^51) . Then projecting out all the states but the first two (|-|-), |— )) we get the following 
two state projected model (TSPM): 



k,{T 
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(31) +Y.^kiflck,a + clj, 



k,{T 

where and are two Pauli matrices in the space spanned by |+) and |— ) and the 
quantities e and A are given by 

(32) e = .g|(+|a + at|-) 

(33) A = ^i{+\a^a\ + {-)-\a^a\-)) 

The latter quantity A is the tunnehng frequency between the two phononic states. A 
similar model has been introduced in Ref. |23| to study the strong couphng limit of 
the Holstein model. The TSPM reproduces exactly the DMFT of the Holstein model 
in two limits: weak coupling and adiabatic limit. In the former case the projection 
of the phonon space has no relevance therefore the TSPM reproduces the perturbation 
expansion developed (in the limit of infinite bandwidth) in Ref. |24j . The adiabatic limit 
is instead recovered as A — > 0. No phonon tunneling occurs and the model can be solved 
exactly by CPA recovering the solution of Ref. |T5j . 

To analytically span from the strong (V^ — *■ 0) to the weak {g 0) coupling 
regimes of the equivalent impurity Hamiltonian it is useful to devise an Iterated Born- 
Oppenheimer Coherent Potential Approximation (IBOCPA) scheme. Starting from the 
Green's function for Vk = H35(l for the spinless case and (|36|l for the spinful one, we 
notice that in both cases Gq(0) can be written ad sum of two contribution Ga(0) = 
(l/2)(Ga.+(0) + Ga,-{0)) where (±) label a phonon state Then we write the prop- 
agator in presence of hybridization as 

(34) G(0) = i y ; — 

where the hopping of the electron from the impurity to the bath is described by the 
cavity-field correlator D^/4G(0) ©. Iteration proceeds substituting the propagator H34|l 
back in Eq. 1261 giving a new BO potential and new TSMP parameters (eq. (I32|l ) and 
finally a new Green's function. Iteration continues after convergence is reached |26| . 

In (|34|l G^^(0) are propagators obtained from the solution of the atomic {Vk = 0) limit 
of the TSPM in both spinless and spinful cases. The spinless atomic Green's function 
can be easily found to be 

a—± ^ 

Ga(w) has a pole at a; = induced by phonon tunneling, whose weight vanishes as A ^ 0, 
accompanied by two resonances at ±21. The zero energy peak is due to transitions in 
which both charge and phonon "spin" change while the side peaks arise from charge 
transfer in a frozen phonon " spin" . In this sense the side peaks are adiabatic features 
which survive when phonon tunneling A — > 0. The spinful atomic Green's function is 



(36) G,(0) = J E ^ 



;-A i + A 



2 ^ 21 \0 + a{l + A) + a{l - A) 
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The most striking difference with, the spinless case (|35|) is the absence of the zero frequency 
pole. In the spinful case the tunneling of the phonon A is always associated to a finite 
energy transition, and it only splits the finite frequency poles associated to the transition 
from singly to the empty or doubly occupied ones. 

IBOCPA assumes that the tunneling states of the phonon in the adiabatic potential 
remain unaltered during an hybridization event. As in standard CPA a band is associated 
to each local level, but our IBOCPA gives a Fermi liquid solution in the spinless case 
for every value of the coupling, as opposed to the case of the Hubbard model. The low- 
energy band arising from the zero energy pole in the zero hybridization limit is indeed 
coherent. This can be easily realized by analysis of the self-energy. When w — > the 
spinless propagator defined by and is dominated by the zero-energy pole of Ga 
(|35|) and consequently the self-energy obtained through (|34|l is purely real. 

Conversely in the spinful case a MIT due to local pair formation occurs in the CPA 
approximation at a critical value of A. Finally we emphasize that we recover the adiabatic 
solution of Ref. ^Hl as 7 ^ for finite A as A ^ 0. In this case the IBOCPA is exact 
and gives the Green's function of Ref. [T3|. However the IBOCPA procedure is certainly 
affected by serious problems approaching the MIT in the spinful case. A more careful 
treatment of the low-energy part of the Green's function has been performed in this case 
in Refs. |23l 1271 IT], where it has been observed a MIT scenario similar to the half-filled 
the Hubbard model, i.e., a quasi-particle peak that shrinks to zero width approaching 
a critical value of A. In the spinless case instead a resonance is present at zero energy 
within a CPA approach. It is not associated to a Kondo effect but rather to phonon 
tunneling which drives charge fluctuations. On the other hand a Kondo like behavior 
can be ascribed to the bipolaron or pair formation. For a discussion of the limitation of 
the IBOCPA approach see also Ref. 25 . 

5'2. Antiadiabatic regime. - While in adiabatic limit the phonon displacement becomes 
a classical variable, and we are left with an electronic model which depends parametrically 
on it, in the opposite limit (7 >> 1) the roles are exchanged, and we have a parametrically 
fixed electronic charge on a given site. In this regime the most reasonable starting point 
is the Lang-Firsov (LF) canonical transformation ^Kl I20| S = exp(r). The generator of 
the transformation reads 



introducing the parameter a — g/ujQ which is the relevant e-ph coupling parameter in 
the anti-adiabatic regime |3UI 15]. 

The canonical transformation diagonalizes the impurity Hamiltonian in the absence 
of hybridization by eliminating the e-ph interaction part. In the spinful case the phonon 
energy term of H14|l gives rise to the well known bipolaronic instantaneous attraction 
[29 | . The hybridization term of H14|) is modified by acquiring an exponential term in the 
phonon coordinates leading to 



(37) 




a 



-T 



(38) 




-2 — (s - l)n^ni V(- - n„) + ujQa^a, 
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where ricr = flfa- Notice that in the anti-adiabatic Hmit 7 — > oo, if A is kept constant a 
vanishes. In this case spinless electrons are not renormalized, while spinful electrons are 
described by an attractive Hubbard model with \U\/D = A. 

If we want to proceed with analytical methods, the hybridization term must be treated 
in an approximate way. Assuming that in the anti-adiabatic limit the impurity density 
is constant during the fast motion of the phonon, we average out the phonon term on 
the displaced phonon ground state. This is the so-called Holstein Lang-Firsov Approxi- 
mation (HLFA), which has not to be confused with the exact canonical transformation 
()37(l . HLFA gives rise to the exponential renormalization of the hybridization constants 
where each Vk is replaced by Vk exp(— a^/2). Such a replacement implies the well known 
exponential renormalization of the bandwidth Dexp(— a^). 

To get the electron Green's function G{lij)^ the explicit action of LF transformation 
have to be taken into account into both creation and destruction operators appearing in 
the definition of the Green function. Following Refs. |H3 and [221, we obtain in both 
spinless and spinful cases 

(39) G{oj) = e-""Gp{u;) + i ^ e-"'^Gp(a; - nujo). 

where Gp{uj) is the Green's function of an impurity (with a negative U interaction in 
the spinful case) with an exponentially reduced hybridization to a bath of conduction 
electrons. The self-consistency condition can be written explicitly in the spinless case 
due to the lack of interaction terms on the impurity: 

(40) G,(^) = -i^, 



4 



G(u;) 



where G(0) is the local Green's function of the lattice. In the spinful case a Lang- 
Firsov Coherent Potential Approximation (LFCPA) can be devised for the resulting 
HLFA attractive Hubbard model with U = — 2g^/cjo giving 

(41) G,(^) = i' ^ ^ 



2\uj~ e-"%G{u;) - C//2 uj - e-"%G{u;) + U/2J' 

Notice that the theory developed here for the Holstein impurity model differs from 
that developed directly in the lattice model[3Tj. In that case an equation identical to H39|) 
is recovered for a band of free electrons therefore giving a low-energy coherent polaronic 
band in the spinless case. It is however easy to show that this form of the spectral function 
is not compatible with a /c-independent self-energy. The self-consistency condition 1)40(1 
gives rise to a non-zero damping at the Fermi level even in the spinless case. However, 
when 7 becomes larger a gets smaller reproducing the anti-adiabatic coherent behavior 
at low energy in the spinless and the negative- 1/ behavior in the spinful cases. 

In the anti-adiabatic regime the HLFA approach gives an estimate of the MIT 

(42) Xmit = \U/D\mit exp(-a2) 



where \U/D\mit ^ 2.94 is the MIT value of the negative-C/ Hubbard model 0. The 
pairing MIT occurs for smaller A as the adiabatic regime is approached (see Fig. EJ. 
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The phonon PDF can be easily derived within LF approach. Being the local electron 
densities parametric variables in the anti-adiabatic we find 

(43) Pix)^J2'^iPo{x-xi) 



where wi is the probability of having an occupancy n; , xi the relative displacements and 
Pq{x) the ground state PDF of an harmonic oscillator. Po{x — xi) is then the conditional 
probability of having a displacement x given a definite occupation n/. In the spinless 
case n; = 0,1 with equal probability giving 

(44) p(x) = l(Fo(x-Xo) + Po(a: + a;o)) 

where xq = ^f\/2. A definite polarization can be associated to the ground state if the 
PDF becomes bimodal. By requiring dP{X) / dX\x=i) > 0, which simply means X = 
turns from a maximum to a minimum, we get the usual anti-adiabatic condition for the 
existence of a polaronic state, i.e., > 1 (see fig. EJ. 
In the spinful case n/ = 0, 1, 2 

(45) P{X) = nd{Po{X - 2Xo) + Po{X + 2Xo)) + (1 - 2nd)Po{X) 

where nd = (ri|n|) is the site double occupancy. It is worth noting that in the insulating 
state Ud — 1/2, and the zero-displacement PDF associated to singly occupied sites is 
depleted. The existence of a definite polarization is now associated to a bipolaronic 
state. The condition under which (|^ becomes bimodal, is 



(46) ea;p(-2a^)(4a^ - 1) > 



2nd 



An estimate for the bipolaronic transition can be obtained by taking nd = 1/2, which 
gives Xpoi = 7/2. The presence of a fraction of singly occupied states increase the critical 
value of A (see Fig. (SJ. We notice that the spinful PDF for = 1/2 maps onto the 
spinless one after the same rescaling of the adiabatic regime H3l)|) . 

5'3. Results from DMFT-ED. - The behavior of P{x) and p obtained from DMFT-ED 
compared with that of the previous theories is shown in Figs. 13 and 01 for adiabatic and 
antiadiabatic regime respectively. In each diagram the value of the quantity shown has 
been shifted upward according to the value of the coupling lambda. The values of 7 and 
A have been chosen according to the scaling (^0)1 . 

Let us first discuss the adiabatic regime. Polaron crossover is seen as a qualitative 
change in the shape of phonon PDF shown in the upper panels of Fig. 13 where BO 
approximation and DMFT results are compared. The anharmonicity due to e-ph in- 
teraction increases as the coupling increases leading first to a non-Gaussian and finally 
to a bimodal PDF at A > Xpoi ■ This behavior signals the appearance of static distor- 
tions, even if we are neglecting any ordering between them. ^From Fig. O is evident 
that BO approximation works well in both the metallic and the polaronic regimes. The 
reason for the accuracy of the BO procedure in the polaronic regime is that, contrary 
to its usual implementation in the weak e-ph coupling |34|. here we take into account 
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the anaharmonicity through Eq. H29() in a non perturbative way. However, BO does 
not accurately reproduce the phonon PDF around the polaron crossover. In this case 
electron and phonon states are strongly entangled, and cannot be approximated properly 
by a disentangled BO state. By a comparison of the spinless and spinful cases in Fig. 0| 
we see that the occurrence of the MIT does not influence much the differences between 
full DMFT and BO, which are in both cases relevant near the polaron crossover. In the 
lower panels of Fig. Owe compare the electronic DOS from ED-DMFT with IBOCPA. 
The different behavior of the spinless and spinful case is not so evident. However, com- 
paring the spinless spectrum with the corresponding spinful, we see that a quasiparticle 
peak is present in the former case, while a depletion of low energy part is much more 
evident in the latter. At strong coupling the discretization of the bath inherent to the 
ED solution of DMFT does not allow us to identify a well defined quasiparticle peak in 
the spinless case. A more careful analysis of the quasi particle spectral weight (B5j shows 
however that no pairing MIT occurs in the spinless case. Notice that the IBOCPA seems 
to be much closer to DMFT-ED in the spinless than in the spinful case where the CPA 
approximation for the electronic degree of freedom is apparently much less adequate. 

The upper panel of Fig. 0] shows that the phonon PDF becomes bimodal at a very 
large value of the coupling. HLFA overestimates polaronicity in the spinless case while 
it behaves better in the spinful case. In the lower panels of Fig. ^the electron DOS is 
compared with the results of HLFA in the antiadiabatic regime. The different behavior 
of spinless and spinful cases is marked here by the presence of a pairing MIT in the 
former, well before the (bi)polaron crossover. HLFA correctly catches the gross behavior 
of the DOS. Notice that at the strong coupling the CPA employed to obtain the lower 
diagrams accurately reproduces both the position and the width of the side bands. 

The different behavior of spinless and spinful system can be easily understood in terms 
of strong coupling anti-adiabatic perturbation theory for the original lattice problem 
introducing the charge sector pseudo-spins 37 . In the spinful case at second order in 
hybridization Vk an anisotropic Kondo Hamiltonian can be obtained (38., .39. 25^ . The 
anisotropic Kondo couplings (see also Eq. (8) of Ref.^Ul) reads 



where the -I- (— ) sign is taken for the |j (_L) coupling and we have assumed for simplicity 
V — J2k^k- In the spinless case the processes leading to J± do not exist while the 
remaining J|| is solely associated to charge fluctuations. As opposed to the charge Kondo 
effect of the spinful case, no Kondo effect is expected for spinless fermions. A strong 
coupling estimates of J|| gives 



and J±/J\^ cx exp(— 20;^) which means an exponential suppression of the superconduc- 
tivity versus charge correlation at strong coupling due to retardation effects |2H1 EHl ■ In 
the anti-adiabatic limit 7 — *■ oo instead the Kondo couplings becomes isotropic. 

The above observations can be summarized in the phase diagrams of Fig. [S] The 
polaron crossover line Xpoi is strongly 7 dependent in both cases. Above this line a 
polaronic (bipolaronic) regime is attained in the spinless (spinful) case. In the spinful 
case, we also have a Xmit line, which separates a normal phase from a paired insulating 




2 



(47) 
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Fig. 3. - DMFT data in the adiabatic regime 7 = 0.1 spinless (panels on the left) and 7 — 0.2 
spinful (panels on the right). The various curves refer to different value of A spanning from 0.1 
to 1.8 in the spinless case and from 0.05 to 1.1 in the spinful case and are shifted according A 
value. Upper panels show the phonon PDF while lower panels the electronic DOS. 




Fig. 4. - DMFT data in the antiadiabatic regime 7 = 2.0 spinless (panels on the left) and 7 — 4.0 
spinful (panels on the right). The various curves refer to different value of A spanning from 0.4 
to 6.5 in the spinless case and from 0.2 to 3.0 in the spinful case and are shifted according A 
value. Upper panels show the phonon PDF while lower panels the electronic DOS. 
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phase [55]. For large phonon frequency we can have pairs without bipolaronic behavior, 
as it can be understood by recaUing that in the antiadiabatic hmit the Holstein model 
becomes an attractive Hubbard model, where no polarization is associated to the pairing. 




- Bipolarons 






Pairs 








Normal 



4 6 8 10 

y 



Fig. 5. - Phase diagrams of the spinless (left) and spinful (right) T = Q Holstein model at half 
filling. Left panel: the bold line is the polaron crossover from bimodality of P{x) and the dotted 
line is the anti-adiabatic estimate \poi = 27 for the polaron crossover. Right panehbold curve 
is the bipolaronic MIT from vanishing of quasi-particle spectral weight Z, thin solid line the 
polaron crossover, bold dotted line is the anti-adiabatic prediction for bipolaronic MIT (Eq. 
(I42ll 'l. light dotted line is the anti-adiabatic estimate A > 7/2 for the polaron crossover. 



6. Half filled Holstein model: spinless vs spinful cases at T > 

Using the QMC procedure described in section I^TI we are able to study the normal 
phase at finite temperature. At fairly high temperature the MIT becomes a crossover 
and we are faced with the problem of finding a suitable quantity to locate unambiguously 
this crossover. In analogy with the phonon PDF used to mark the polaron crossover we 
can define a distribution of a quantity that locates the pairing crossover. Let us define 
the distribution of the center of mass Xc ( "centroid" ) of the phonon path 




where the averages are evaluated over the action In the same formalism the 

phonon PDF is the distribution of the the endpoint x = a;(0) — x{(3). 

The meaning of the centroid variable Xc has been discussed in for a single 
particle in a binding potential. Here the variable X represents the fluctuating position 
of the particle, and Xc is the classical position of the particle For an heavy particle, 
the classical limit holds, so that P{X) and P(Xc) coincide The lighter is the 

particle, the broader the wave function, increasing the variance of P{X) while P{Xc) 
turns out to be essentially determined by the binding range of the potential. Here, we 
use P{Xc) for the many-body problem, and propose that pairing can be associated with 
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a multimodal behavior in P{Xc) which takes place at a given value of the coupling Xpair- 
(¥2) The ability of our estimator to determine the pairing crossover can be understood 
by inspecting the interaction term 1)13(1 . In the adiabatic limit (7 — > 0) the kinetic term 
forces the phonon path to be r-independent and the phonon field becomes classical. In 
this limit Xc is equal to X and P{X) and P{Xc) obviously coincide. Thus the centroid 
distribution becomes bimodal when the system is polarized, i.e., Xpair — Xpoi, which 
is exactly what one expects since a static field can induce pairing only with a finite 
polarization. Notice that in this sense the bimodality of P{Xc) is a precursor of the 
actual pairing MIT which occurs at T = and o^o = at a larger value of the coupling 
[El (see fig. EJ. 

On the other hand, in the opposite atomic (Z) — + 0, 7 — + 00) limit the electron 
density becomes a constant of motion. Therefore Eq. (|13|l takes the transparent form 
— \/A(n — 1) cItx^t) where the electron density is directly coupled to the centroid Xc- 
The average appearing in Eq. H49() is readily carried out, giving 



(50) P{Xc) oc exp 



- ^ log(2cosh(/3VAXe + 1))) 
2 fi 



which becomes bimodal for A > Xpair — ST. This is exactly the scale where double 
occupancies start to proliferate in the atomic limit. Therefore the bimodality of P{Xc) 
correctly signals the onset of pairing also in the antiadiabatic regime. In the same limit, 
it can be proved that the endpoint distribution P{X) has a variance which scales with 
1/VAt and as a consequence no definite polarization may occur. We finally notice that 
the D ^ limit of adiabatic P{x) 15. coincides with P{Xc) of Eq. (|50fl . Since for uq = 
the distributions of X and Xc coincide, we conclude that in the atomic limit P{Xc) is 
the same for loq — and loq — 00. This suggests that the pairing crossover may depend 
on Wo more weakly than the polarization one. 

6-1. Results from DMFT-QMC. - To analyze the evolution of P{X) and P{Xc) at 
finite D and wo we use DMFT-QMC. The numerically exact results, shown in FigEl 
clearly show that P{X) and P{Xc) tend to coincide in the relatively adiabatic case 
7 = 0.1, as expected from the previous arguments about the adiabatic limit. The two 
quantities are clearly different for 7 = 1 and 8. For temperatures smaller than ojo, the 
polarization crossover Xpoi moves to larger values as 7 is increased, while the line (Xpair) 
where P{Xc) becomes bimodal is only slightly shifted to larger couplings with increasing 
7. This is strongly reminiscent of the behavior of the metal- insulator transition in the 
Holstein model at T = 0, whose critical coupling is slowly increasing with 7 and then 
saturates to the asymptotic 7 = 00 value (SS). The polarization crossover is instead 
roughly proportional to 7[3S|. Both at zero and at finite temperature the line where the 
centroid becomes bimodal does not coincide with the metal-insulator line, but it can be 
considered a precursor which depends in a very similar way on 7 and on A. 

Our DMFT results can also be compared with the semi-analytical results for 7 = 6 
[T^ (thin dashed line in FigEJ. The 7 = 0.1 case is in very good agreement with the 
adiabatic result, and also the 7 = 1 and 8 cases, at high temperature, fall on the same 
curve. The centroid distribution depends weakly on 7, as suggested by the atomic limit. 
Interestingly, the adiabatic result displays a re-entrance at low temperatures as a non 
monotonic behavior of Xpair and Xpoi as a function of temperature. Although the QMC 
simulations do not reach sufficiently low temperatures, we find that the re-entrance is 
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Fig. 6. - \pair and \poi at 7 = 8. The dashed hne represents the 7 = result where Xpair = Xpoi- 
The dcished arrow indicates the zero-temperature results for the polaron crossover Apoi for 7 = 8. 



present also for 7 different from zero, as indicated by the arrows in FiglBJ which mark 
T = results for the polarization crossover for the Holstcin model [53]. 

7. Conclusions 

The normal state properties of strongly coupled electron phonon systems can be drawn 
from the comparison of electronic spectral and phonon PDF properties. 

Qualitative changes in phonon PDF signals a polaronic crossover while an electronic 
MIT can be seen from a gap in the electronic DOS or by the vanishing of the quasi 
particle spectral weight. This last transition can be observed only at sufficiently high 
density provided the Coulomb repulsion is neglected. In the limit of low density the 
carrier show a polaronic behavior trough developing a definite polarization in the phonon 
PDF while at large densities occupied and empty sites makes the phonon PDF bimodal. 
At the polaron crossover fluctuation of phonon coordinates tends to be larger than those 
at any other coupling. In this regime Born-Oppenheimer approximation is shown to fail 
also when phonon frequency is much less then electron bandwidth. For a single hole in 
the t— J model a further source of localization is due to magnetic superexchange energy 
which tends to localize the spin-defect spreading due to the presence of the hole. 

At non zero temperature the pairing MIT becomes a crossover. To locate the point 
of this crossover it is possible to define a suitable distribution associated to the phonon 
classical position which is coupled to the electron density in the non-adiabatic regime. 

There are several limitations in the DMFT method used to obtain the present results. 
The biggest one is the non dispersive nature of bosonic excitation. In the t — J case 
the absence of spin-wave dispersion lead to localization of the single hole. However local 
properties seem to be quite well represented by DMFT in when compared with those 
found at finite dimensions [H]. In the half filled Holstein model case bipolarons do not 
have coherent motion therefore the vanishing of the single particle spectral weight implies 
a MIT. Several extensions of DMFT such as Extended-DMFT for the < - J model HH] 
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or CDMFT for the Hubbard model try to overcome such hmitations. Work along 
these lines for the Holstein and Holstein-Hubbard model are currently in progress. 
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